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Abstract - We derive the stochastic equations of motion for a tracer that is tightly attached to a 
semiflexible polymer and confined or agitated by an externally controlled potential. The generalised 
Langevin equation, the power spectrum, and the mean-square displacement for the tracer dynamics 
are explicitly constructed from the microscopic equations of motion for a weakly bending wormlike 
chain by a systematic coarse-graining procedure. Our accurate analytical expressions should 
provide a convenient starting point for further theoretical developments and for the analysis of 
various single-molecule experiments and of protein shape fluctuations. 
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Introduction. — Many of the tools available to an 
experimental biophysicist probe either the fluctuations of 
I semiflexible polymers or their response to external forces. 
'^This includes dynamic light scattering (l]|2], active and 
S passive microrheology of polymer networks^ or cells jl] , 
^magnetic bead twisting cytometry [s], DNA relaxation 
I— &nd stretching [6], single-molecule force spectroscopy [t] 
and electron transfer techniques (s] . From a theoretical 
point of view, these methods expose different aspects of 
the wormlike chain (WLC) model. It provides a minimal 
l/^ description of semiflexible polymer physics in terms of 
^/~) an inextensible, thermally fluctuating elastic beam, and 
'^^has found broad acceptance on grounds of its excellent 
00 agreement with experimental data. 

o 

Deformations of the polymer contour excite a broad spec- 
. . trum of bending modes with widely disparate relaxation 
^ times, thus resulting in anomalous, subdiffusive dynamics. 
In the linear regime, valid for equilibrium fluctuations or 
small external forces transverse to the polymer backbone, 
C^the dynamic mean-square displacement of a tagged (but 
mechanically unaltered) monomer obeys MSDj_(i) cx t^/'* 
9|[T0 . For fluctuations parallel to the polymer axis, the ad- 
ditional longitudinal solvent friction induces tension forces 
which in turn stiffen the polymer and give rise to a dif- 
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ferent scaling behaviour, MSD||(t) oc t"^/^ 
also dominates the response to strong point 
externally imposed solvent flows; the resulting equations of 
motion are highly nonlinear and can produce a multitude 
of different dynamical regimes even in the course of a single 
experiment 



In many cases of practical interest, a full evaluation of 
the dynamics would be needlessly complicated and waste- 
ful, since experimental manipulation and data acquisition 
are strongly localised, say, to an attached tracer particle, 
or a tagged monomer, in the following simply referred to 
as "the tracer". Farther parts of the polymer matter only 
insofar as they contribute to the force on the tracer. It can 
then be preferable to integrate out the polymeric degrees 
of freedom beforehand and subsume them under an effec- 
tive equation of motion describing the tracer coordinate 
only. Such a reduced description is for example known 
for the important special case of a tracer subjected to an 
externally prescribed deterministic force protocol 16 . It 



cannot, however, easily be extended to accommodate for 
the ffuctuating forces exerted onto the tracer by an exter- 
nally controlled confinement potential. Practical examples 
that involve such a potential are provided by various single- 
molecule manipulation techniques (think e.g. of an actin 
filament labelled with a gold nanoparticle that is trapped 
by optical tweezers). The analysis of high-frequency shape 
fluctuations of globular proteins, as measured by electron 



transfer techniques |8 17 provides another important ex- 
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ample. Indeed, the WLC has been proposed as one possible 
model of protein fluctuations, but to date only numerical 
evaluations of the corresponding noise and friction func- 
tions are available within a mean-fleld approximation to 
the WLC (20]. 

In the following, we systematically derive the sought- 
after reduced equation of motion for the tracer coordinate 
x{t), starting from the WLC in the weakly-bending limit. 
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which is asymptotically exact for large bending rigidity or 
strong stretching force. The resulting GLE, 



[ ATK(t-T)x{T)=F{. 

Jo 



x,t) + E{t) 



(1) 



is not necessarily linear, as it may include an arbitrary 
external potential U{x, t) = ~ J dx'F{x' , t). In the infinite- 
length limit L — oo, the GLE ([ij can be worked out explic- 
itly in terms of the microscopic parameters, which comprise 
the length L of the polymer, its bending rigidity k, and 
tension /. It is validated by comparison with numerical so- 
lutions of the exact equations of motion. We moreover give 
a simple interpolation formula ([9| that provides a universal 
description for polymer-bound tracer particles in strongly 
localised externally controlled potentials and possibly also 
for the mentioned protein shape fluctuations. We expect 
it to become a valuable and convenient tool for analysing 
a wealth of experimental data and for future theoretical 
developments. By providing a physically transparent and 
concise parametrisation of measured tracer movements, 
it will moreover be helpful in the mutual comparison of 
data obtained with diverse experimental techniques. To 
exemplify our approach, we calculate various observables 
characterising the time-dependent spatial correlations of 
the tracer motion, such as its power spectrum and mean- 
square displacement in presence of a harmonic trap. 

Langevin description of a WLC. In the WLC 

model, a semiflexible polymer is mechanically represented 
as an inextensible elastic beam of length L and bending 
rigidity k. It follows that thermal forces can only induce 
significant bending on length scales larger than the per- 
sistence length £p — K,/[k^T) (in 3 dimensions). In the 
weakly-bending limit, valid for large persistence length 
ip ^ L 01 strong external stretching force / ^ k^T/ip, 
the polymer is essentially straight and can thus be treated 
in terms of its small excursions r±{s, t) from the straight- 
rod ground state. The elastic bending energy in a given 
configuration r±{s,t) reads |9j 



n 



ds 



Shape fluctuations of the polymer then obey a Langevin 
equation obtained by balancing the corresponding bending 
forces with friction and thermal (Gaussian white) noise |9], 
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^~nr'^' + frl+U (2a) 
= (2b) 
{U{s,t)U{-'^'t'))^2(^^kBT5{t~t')5{s~s'). (2c) 

In the vein of a similar treatment for a Rouse chain 




Figure 1: Force diagram for the combined system of a tracer 
(red) and an attached WLC (blue). The tracer is displaced by 
an external potential U{x,t) = — J F{x, t) Ax. The constrain- 
ing force pair ±7(t) fixes the polymer backbone to the tracer 
position X. 



it absorb the external driving force F{x,t), 

Ctri(t) =F(x,t)-7(t)+etr(i) (3a) 

+ 7(i)<5(s-so)+e±(5,t) (3b) 
r±{so,t) = x{t). (3c) 

Here Ctr and ^tr denote an optional friction coefficient and 
Gaussian white noise source for the tracer, respectively. 
Both may be set equal to zero in the tagged- monomer 
case. The (Lagrange) forces ±7(t) serve to rigidly tie the 
tracer to the polymer contour at Sq, as required by the 
constraint ( 3c 1. The latter is simplified by transferring to a 



comoving reference frame, in which the tracer is always at 
rest, A(s, t) = r±(s, t) — x{t). The corresponding equation 
of motion acquires a spatially constant friction term due 
to the drift between the comoving reference frame and the 
solvent, 

aA(s, t) = -A^A""(s, t) + fA"{s, t) + e±(s, t) 

+ l{t)5{s - so) - Ci^x{t) (4a) 
A(so,t) = 0. (4b) 

The formal solution for 7 can immediately be written down 
in the frequency domain. 



/ dg ^_L,c^((T)G',^(so,g) 
Gu{so, So) 



/ do- ^.^(so,^) 
Go; (so, So) 



monomer 21 , we now introduce a tracer at s = sq and let 



where G{s,s',t) denotes the transverse Green's function 
of a WLC, i.e. its transverse deformation r±{s,t) in re- 
sponse to a unit force impulse S{s — s')S{t). This explicitly 
establishes the link between the microscopic equations of 
motion ^ and the coarse-grained equation for the tracer 
([ij. A practical way of evaluating 7(i) numerically consists 
in the decomposition of A{s,t) into its normal coordinates 
|22j. This procedure is discussed for arbitrary boundary 
conditions and monomer positions sq in the appendix. 

Analytical solution. — To proceed analytically, we 
now consider the centre monomer only, so = 0. Instead 
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of explicitly including the "adhesion force" j{t)S{s), which 
induces a coupling of different eigenmodes and thus renders 
the dynamics nondiagonal, we can then easily take care of 
the constraint A(0,t) = by restricting the function space 
accordingly. Using only those eigenmodes Wn{s) satisfying 
yV„(0) = 0; the singular force 7 can be read off as follows. 

K{A"'iO+,t)-A"'iO-,t))=jit). 

Since the above expression vanishes for odd eigenmodes, 
only even modes Wn(s) — W„(— s) contribute to 7. Re- 
quiring further that A"" constitutes the highest (and only) 
singularity, we thus find A'(O^) = and so obtain the 
following friction kernel. 



L/2 



W„(s) ds 



where 

«wr(s)-/w;'(s) = £„w„(s) 

w„(o) = w;(o) = o, 

and the outer boundary conditions are dictated by the 
physical situation. In the infinite-length limit L — >■ 00, 
valid for i <C ti, K{t) becomes independent both of sq and 
of the choice of outer boundary conditions. Using torqued 
ends yV4(L/2) = yy,7(L/2) = for convenience, we find 



K{t) - 8 



00 

E 

n=0 



L3 



-t/r„ 



(5) 



with a relaxation time spectrum 



K7r4 (2n)4 + (2n + l)2///L- 

Here, /l = ktt^/L^ denotes the Euler buckling force. For 
semiflexible polymers, /l is usually small in comparison to 
externally applied stretching forces. Mode numbers below 
a critical value ric = (///l)^^^ are then tension-dominated, 
whereas shorter-wavelength modes exhibit force-free relax- 
ation and therefore t„ oc ti"**. This divides the frictional 
response K{t) to a transverse force into three different 
asymptotic regimes. At short times t ^ = C-Lk/2/^, 
the backbone tension can be neglected so that K{t) sim- 
pUfies to 

2^2^1/4^^ ^ 



(6) 



r(i/4) 

Between and the terminal relaxation time ti, we have 

2 



r(l/2) 

and at very long times i ^ ti , the exponential relaxation 



K{t > Ti) cx 



(8) 




Figure 2: The analytical interpolation formula ||9| compared to 
the exact numerical solutions for / = (O) a-nd / = IO^/l (□). 



of the lowest mode provides a natural physical cutoff of 
the scale-free intermediate asymptotic dynamics. Using 
a time-dependent mode cutoff at r„ = t, we arrive at the 
approximate interpolation formula 



K{t) 



2V2 



(9) 



which faithfully reproduces the general solution described 
in the appendix, see fig.([2|. 



An exact treatment of the binding constraint ( 3c 



malises the terminal relaxation time, which is why we 
consider t* as a free (fit) parameter. In the long-polymer 
limit (L,r* 00, e^*/'^* 1), eq. ^ reduces to a two- 
parameter formula. 

The associated colored noise term S(i) of eq. ([T]) can 
be determined from the fluctuation-dissipation theorem 
(FDT), 

{Em{t'))=2kBTK{\t-t'\). 

Example applications. — As an example application 
of our single-coordinate equation of motion eq. (jlj , we first 
rederive the time-dependent MSD of a monomer in the 
bending-dominated regime. We then include an external 
harmonic potential to compute both the time-dependent 
MSD and the power spectrum of a polymer-bound tracer 
particle held in a harmonic trap. 

For the free polymer in solution, we set F{x,t) = and 
hold the tagged monomer fixed until t = 0. Using the 
early-time asymptote to K{t), eq. ([6|, its trajectory then 
follows as 



(7) x{t) - x(0) 



1 



* S(T)dT 

2^r(3/4)Ki/4^f 4 (t-T)i/4' 



(10) 



The MSD follows from eq. m 



MSDF=o(i<r„J 



1 V^fceT 
2r(7/4) k1/4 



t 



3/4 
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oj [Hz] 

Figure 3: Top panel: power spectrum of an unmodified tagged 
monomer (no excess friction) with ip = 10 /im, = lOmPa ■ s, 
k = IpN/nm, T = 300K, L ~> oo {i.e. ujt, oo) and / = 5pN 
(□), / = (O)- Bottom panel: power spectrum of an at- 
taclied tracer particle that has a perceptible friction coefficient 
(^tr = BTrr^waterJ", r = 10 um, other parameters as above. The 
Lorentzian power spectrum of the same bead without the at- 
tached polymer is shown for comparison (dashed line). 



This subdifFusive behaviour coincides with previous the- 
oretical predictions [2] and has been measured directly 
and indirectly in networks of polymerized actin [T 10 and 
microtubuli [5] . Including a stationary harmonic trap of 
stiffness k, i.e., F{x,t) = F(x) =^ —kx, equation (10 1 turns 
into 



x{t) - a;(0)E3 



where 



E„(2) 



^ r(i 

n=0 ^ 



an) 



(11) 



denotes the Mittag-Leffler function, which can be regarded 
as a generalised exponential function. The resulting MSD 
for x{0) = reads 



MSDF^o(i<T„J 



kjiT 



1 - El 



272^1/4 



A deterministic external force fextit) or a dynamically 
moving optical trap, represented by /ext(i) = kxo{t), can 
be included along the same lines by setting F{x,t) = 



—kx + fcKt{t). Finally, we calculate the power spectrum 
density for a polymer-bound bead in harmonic confinement, 



{x{t)x{0))e-''^'dt , 



as measured by video microscopy or other in-plane tech- 
niques. (Note that for 3c? position tracking x becomes 
a 2d vector transverse to the polymer backbone, and 
^uj^ ~ 2S'tj.) 

The somewhat lengthy explicit result has a simple struc- 
ture, suitable for fitting experimental data. Specialising to 
the infinite-length limit {ujt^ — oo), we find 



k^ 



(12) 



where R^^ and /^^ denote the real and imaginary parts (with 
a dimension of length) of 

2V2^[{J - ^ p + 4iLuC±K)'^^^ 

+ (/+V<F+4^^)-i/2], 

respectively. For a tracer that causes a perceptible friction 
(Ctr = 67r?7rtr > 0), the ultimate high-frequency limit is of 
the usual Lorentzian form S{uj — > oo) ^ uj^^. If (tr van- 
ishes or is imperceptibly small compared to the monomer 
friction, the decay of the power spectrum is slightly weaker 
at large frequencies, S ^ uj^'^ 1^. This can intuitively be 
understood as a consequence of the frequency-dependent 
"apparent bead size", given by the subsection (of length 
where ~ C,±^(jj/k) of the attached polymer that 
equilibrates with the bead within one period uj~^. See 
fig. [3] for a graphical representation. 

In the absence of backbone tension (/ = 0), the power 
spectrum (12) simplifies to (s — sin(7r/8), c — cos(7r/8)) 



sL 



=0 



8Cf"(«MV4(,s + e) + 2Ct. 
[fc + 2Ki/4(ac^)3/4(s_c)]' 

+ uj^ [Ctr + 2Ci/'(^Mi/4(s + c)' 



Conclusions. — Starting from the formally exact 
WLC equation of motion in the weakly-bending rod approx- 
imation, we have derived a generalised Langevin equation 
describing the dynamics of a tagged monomer of (or "tracer" 
attached to) a semiflexible polymer. The tracer was al- 
lowed to be under the influence of an arbitrary external 
potential. Our method is simple, direct and analytically 
solvable. We have furthermore derived a uniformly valid 
analytic interpolation formula which may serve as a com- 
pact (two- or three-parameter) parametrisation of, e.g., the 
motion of a tracer attached to a semiflexible polymer and 
manipulated by an optical trap, or of conformational fluc- 
tuations of protein domains. With regard to quantitative 
applications using metallic tracer particles in combination 
with optical traps, it might be worthwhile to extend our 
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results along the lines of f23l to take the heating of the 
tracer into account. 
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Appendix. — For the evaluation of j{t), we decom- 
pose A(s,t) into bending eigenmodes Wnis) and mode 
amplitudes an{t), such that 



N 



where Afi_2 and Vx^(^ are given by 



W2(so) -Wi(so) 



Wiv(so) -Wn_i(so) 
W„(so) 



Wo(so) Wi(so) 
Wi(°o) Wo(°o) 

_W2l£0) Wl(°0) 



Wjv(so) Wjv_i(3o) 





A(s,t) = ^W„(s)a„(i) 



K = 



nW:'\s)-fW:{s) 

L 

daW„(a)W„,(a) 



5n 



The actual physical boundary conditions at the polymer 
ends dictate the detailed functional form of the yV„. A 
free polymer in solution requires |24| 

w::(-L/2) = w:'{-L/2) = w:{Lm = w;"(l/2) = o, 

but the precise choice of boundary conditions is irrelevant 



AiWo(so)- Wi(so)^o 



AnWn~i{so) - Wn{so)An-i 


eoWWi(so)-Wo(so)eiW 



1 



^N-i{t)WN{so) " Ww-i(so)CwW 




The solution to eq. ( 13 I then reads 



to the following discussion. Projecting eq. (4b I onto each 



of the Wn , we find + 1 distinct equations of motion for 
the a„ , 



+ rdre^^i"^^^(*--)Mri(i(T)14 + V^ir)) . 

J to 

Inserting the above solution into the original equation of 



d„(t) = -a„{t) — + i{t) (W„(s) I 1 

{Wn{s) \ Ui-S,t)) , Jit) 



Ci 



+ ^W„(so) ) , (A.l) 



where t„ is the relaxation time of the nth eigenmode. 
The singular force term can be eliminated by choosing 
a complete set of allowed displacements in mode space, 
I]>V„(so)a„(i) = 0, leading to 

51 = (Wi(so),-Wo(so),0,...,0) 

52 - (0,W2(so),-Wi(so),0,...,0) 



motion (A.l), we obtain the force of constraint in mode 
space, 

= C±dtan{t) + (±anit)/Tn - £,n{t) + i{t)A.nC± 



/ dr{Mr'M,+dmg{r-'))e''^ *^^(*-) 

Jta 



5n = (0,...,0,Wjv(so),-Ww-i(so)) 



(A.2) 



(A.3) 



The part outside the integral vanishes, except for the A^th 
entry, which is an artefact introduced by the mode cutoff 
and can safely be ignored: having a finite minimum bending 
mode wavelength implies that even at arbitrarily short 
times, a finite part of the polymer will be dragged along 
with the monomer. We thus obtain ^{t) as follows. 



Combining the eqs. (A.l I of motion with the corresponding 



constraint equations (A.2), the dynamics of the a„ are 



determined completely, reading 

Midt{ao, . . . , anV = A'hiao, . . . , unV + x{t)Vx + V^it), 



lit) 



lim 



n=0 -I 



+ diag(r,7i))e*-^^"^^^(*-^)Mr^(i(T)V; + ^(r)) , 
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where W = C-l(Wo(so), . . . , Wjv(so))- Identifying the ran- 
dom and the i-dependent terms with S^{t) and K{t), re- 
spectively, we find 

K{t) = W^{M^^M2 + diag(r-i))e*^i"^^*Mf V, , 
m = f dt' {M^'M2 + diag(T-i)) 

J to 

Note that if sq coincides with the center of the polymer, 
antisymmetric modes will not contribute to "f{t); the calcu- 
lation then has to be restricted to the symmetric component 
As{s,t), otherwise Mi would be degenerate. The above 
procedure trivially extends to an inhomogcncous stationary 
force profile / = /(s). In that case, the differential oper- 
ator ndg — fdg turns into ndf — f'{s)ds — /(s)i9|, which 
changes both the eigenmodes yV„(s) and their respective 
eigenvalues, but the calculation in mode space remains 
unaffected. 
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